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ABSTRACT 

We present a direct image of the innermost companion to the red giant o Andromedae 
using the Stellar Double Coronagraph at the Palomar Observatory. We use a Markov- 
chain Monte Carlo based algorithm to simultaneously reduce the data and perform 
astrometry and photometry of the companion. We determine that the companion is 
most likely a main-sequence K-type star and is certainly not the previously hypothesized 
white dwarf. 


Introduction 


6 Andromedae (K3 III) is red giant with a visual magnitude of 3.28. It has a UV excess which 
implies a hot, high-velocity wind and a 60 and 100 qm excess (Judge et al. 1987) which is most 
likely due to a debris disk (Decin et al. 2003). It is the brightest star in a quadruple system; of the 
outer companions (28.7 and 48 arc seconds), the first has been classified as an M2 V star with V = 
11.3, probably physically associated with the primary as it shares the same proper motion (Bakos 


1976). The outer component does not share the proper motion of the system and is most likely a 


background object (ibid). 


6 Andromedae is a spectroscopic binary with a rather long period of about 57 years (Massarotti 


et al. 2008); see Table [I] for a summary of its physical properties. The presence of the secondary 
has been confirmed photocentrically and astrometrically (Gontcharov & Kiyaeva 2002), with both 
sources deriving an eccentricity of about 0.5. The companion has been conjectured to be either a 


main-sequence star later than G-type (Judge et al. 1987) or white dwarf near the Chandrasekhar 


limit (Gontcharov & Kiyaeva 2002). It has never been directly imaged, however, due to the sec¬ 


ondary’s faintness and proximity to the primary. In this work we image the companion for the 
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first time and measure its magnitude and separation. We find that the secondary is at approxi¬ 
mately the expected separation and determine that it is not a white dwarf, but most likely a main 
sequence star of K-type. This work demonstrates the potential of high contrast imaging with low 
inner working angles applied to spectroscopic binaries. 

Physical properties of 5 And 


Mass (M^+Mb) 

2.6 ± 0.4 M 0 

Gontcharov & Kiyaeva 

( 

2002 

) 

A4+&B 

0.62 ± 0.04” 

Gontcharov & Kiyaeva 

( 

2002 

) 

Period 

57.6 ± 1.09 yr 

Massarotti et al. 

( 

2008 


Radius 

13.6 ± 0.3 R 0 

Piau et al. 

0 

>011 

) 



Luminosity 

68 ± 4 L 0 

Piau et al. 

0 

>011 

) 



Surface gravity (log g) 

2.0 ± 0.3 

Judge et a 


(1987 

) 


Temperature 

4315 ±9K 

Massarotti et al. 

( 

2008 

) 

Metallicity [Fe/H] 

0.04 dex 

Massarotti et al. 

( 

2008 

) 

Rotational velocity (y sin i) 

6.5 km/s 

Massarotti et al. 

( 

2008 

) 

Age 

3.2 Gyr 

Decin et al. 

(2003 

(v- 

uncertain) 

Parallax 

0.032 ± 0.001” 

van Leeuwen 

(200 

7) 



6 And b properties (this work) 
AM (Bracket-y) 

Angular Separation 
Position angle 
Physical Separation 
Spectral type 


6.22 ± 0.05 
0.357 ± 0.0035 ” 
56 ± 1° 

11.55 ± 0.13 AU 
K4 ± 2 


A c = 2.18 n m, AA =0.03 /mi 


parallax from Hipparcos , as above 
derived from ATLAS9 spectra 


Table 1: Previously measured properties of 8 And and newly measured properties of the companion 


2. Instrumentation, Observations and Data Analysis 
2.1. Instrumentation 

The Stellar Double Coronagraph (SDC) is a JPL-developed instrument designed for high- 
contrast imaging of close-in companions to stars, particularly exoplanets. It uses two optical vortices 
in series to simultaneously diffract starlight out of the pupil of the instrument and partially correct 
for the secondary obscuration of the telescope (Mawet et al. 2011b). It has an inner working angle 
of approximately 1A/D, or 90 mas in K-band (2.2 jam) behind the 5 m Hale telescope. It is installed 
between the P3K adaptive optics system (Dekany et al. 2013) and the near-IR imager PHARO 
( |Hayward et al. 2001). 











































































2.2. Observations 


We observed 8 And on October 8-9 2014 (UTC), during the course of normal science operations. 
The seeing was 1.2”, with the adaptive optics system delivering a Strehl ratio of about 85% at an 
airmass of 1.02-1.03. Our observing strategy involved frequently dithering between the target star 
and a reference star, then using post-processing to subtract the speckle pattern from the target 
images using the reference (see the following section for more detail). A Bracket-y filter was used 
concurrently with neutral density filters to reduce the flux from the target when off the coronagraph 
to below detector saturation. Absolute transmissivity of the neutral density filters was measured 


separately, and found to be consistent with Metchev & Hillenbrand (2004). Sky backgrounds were 
interspersed with the target and reference star observations; sky flats were taken five days later. A 


summary of the observations is presented in Table 2.2 


Observing date: Oct. 9-10 2014, JD 2456939-10 


Target 

Images 

Filters 

Exposure Tinre[s] 

Purpose 

6 And 

29 

Br-y, ND2 

9.91 

Photometry, Astrometry 

(3 And 

100 

Br-y, ND2 

2.83 

Photometry 

8 And 

10 

Br-y, ND3 

2.83 

Non-coronagraphic, Photometry 


Table 2: Summary of observations. 


2.2.1. Data Analysis 


One of the main challenges in high contrast imaging is trying to remove speckles due to 
aberrations in the optics after the wavefront sensor. There are a number of ways to tackle this 
contrast-limiting/quasi-static aberration problem; our strategy is sometimes called “reference dif¬ 


ferential imaging” (Mawet et al. 2011a), which involves dithering between the target and a nearby 


star of similar visible magnitude, spectral type, and airmass. This leads to a similar AO correction 
and gravity vector, ensuring a similar speckle pattern. It is then possible to remove some of the 
speckles by either subtracting the reference image or using a more advanced image processing tech¬ 


nique such as the Karhunen-Loeve eigenimage decomposition (Soummer et al. 2012). The latter 


method gives better results than the former in terms of contrast, but has the unfortunate side ef¬ 
fect of reducing the flux of any nearby companions that might be in the image, therefore rendering 
accurate photometry difficult. In this paper, we use a slightly different approach where we forward 
model the target image as a combination of a scaled reference image and a shifted, attenuated 
point-spread function image. This method has some advantages that will be explained below. 


We acquired coronagraphic images of 6 Andromedae and the reference star, (3 Andromedae. 
We aligned and median combined these images after flat-fielding, background subtraction and bad 
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pixel removal. We derive a relative magnitude and offset between the star and companion using 
a Markov Chain Monte Carlo (MCMC) fitting algorithm (Foreman-Mackey et al. 2013). This is 


somewhat different than the usual approach to analyzing fluxes and positions, where one prioritizes 
maximising the signal to noise ratio of the companion, often performing astrometry and photometry 
separately. Here the image reduction, raw photometry, and astrometry are all performed at the 
same time by the MCMC algorithm. There are a few advantages to doing everything at once with 
MCMC. First, one can measure the precision of the reduction algorithm much more accurately: 
the per-pixel uncertainties are Poissonian and straightforward to propagate in the model above. 
Furthermore, the MCMC returns marginal likelihoods, which shows the precision in each parameter 
as well as any correlations. Finally, one does not need an analytic model of the PSF but can use 
images of the instrumental point spread function taken off the coronagraph. This reduces the 
number of parameters in the model, decreases degeneracy, and improves accuracy. 


The generative model for the image data is 


T[x, y] = R a ■ R[x , y\ + P a - P[x -x c ,y- y c \ 


(1) 


where T is the coronagraphic image of 6 And, R a is a constant scale factor, R is the coronagraphic 
reference image of f3 And, P a is another constant scale factor, and P is the point-spread function 
(ie, a unit-intensity normalized, non-coronagraphic image of a point source). Images T, R, and P 
are all single median images. The indices x, y refer to pixel coordinates, and the factors x c , y c are 
shifts in point-spread function imaging data (ie, P[x-1, y-0.34] corresponds to a pixel shift of 1, 
0.34). The constant R a is to correct for the fact that the background speckle field in the science 
image is of a different mean intensity, due to differing stellar magnitudes. P a is the intensity scaling 
prefactor of the point spread function of the companion. MCMC is used to solve for x c , y c , P a , 
and R a simultaneously; the results are presented in Figure [lj The “reference subtracted” image, 
T[x,y] — R a ■ R[x,y] is shown in Figure [2j 


In order to determine the relative brightness, we similarly use the unit-intensity PSF model to 
fit a non-coronagraphic image of 6 And, and the derived intensity allows us to establish a relative 
intensity in magnitudes. The uncertainty in relative intensity is dominated by the uncertainty on 
the neutral density filters’ absolute extinctions. For the companion location, the typical error in 
this case for x c and y c was about 0.01 pixels, or less than a milliarcsecond at 0.025”/pixel. However, 
this is not the true uncertainty in separation because the primary star’s image is suppressed and 
distorted by the coronagraph and its true position is not obvious to calculate. In order to locate 
the position of the primary, we imposed a waffle pattern on the deformable mirror of the adaptive 
optics during observations. The waffle generates astrometric spots 3.9” away from the primary, 
which can be used to locate the position of the star, and we verified our result using the Radon 
transform (Pueyo et al. 2014). The waffle centration has an uncertainty of about 0.1 pixels, which 
dominates the total separation uncertainty. 
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Fig. 1.— Left: a) the background-subtracted target median image, b) the background-subtracted 
reference star median image, c) the background-subtracted point-spread function image, d) best¬ 
fitting model from the MCMC algorithm combining images b) and c) attempting to match a) as 
explained above. The stretch is nonlinear to better show the companion and speckles. Right: 
All the one and two dimensional projections of the posterior probability distributions of the pixel 
shifts (x c , y C: the reference background scaling factor ( R a , and the PSF amplitude used to fit 
the companion P a . The two-dimensional projections show very little covariance among any two 
parameters, and the marginal distribution histograms (along the diagonal) are nicely peaked. 
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Fig. 2.— The reduced, background-removed coronagraphic image of 8 Andromedae. The first Airy 
ring is visible around the companion. The stretch in the image is linear. The colorbar shows the 
relative intensity (as a fraction) compared to the primary 
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3. Results and Conclusions 

The results of the above analysis are shown in Table [l] Judge et al. suggest that the secondary 
companion is either a main sequence star later than G type or a white dwarf. Gontcharov and 
Kiyaeva measure a mass fraction ms/ijnA + mB) = 0.5 ± 0.1 for the binary system and then favor 
the white dwarf assumption, suggesting an even split of mass between primary and secondary, 
placing the white dwarf very near the Chandrasekhar mass limit. 



Fig. 3.— Comparison of the approximate fluxes of 8 And A, three K dwarfs, and a 50000K white 
dwarf. The spectral models are from Castelli Sz Kurucz (2004). The width of the shaded bar is the 
span of the Bracket-y filter. 


However, assuming that the companion is a white dwarf, its radius is constrained to be about 
0.01 Rq, as white dwarfs of 0.5-1.4 M 0 span the radius range of 0.014 to 0.005 R 0 . Comparing 
the expected flux levels of a hot white dwarf to that of the primary (see Figure [3]), the magnitude 
difference through the Bracket-y filter would be about 12 magnitudes, not the measured 6, a 
discrepancy of greater than 100 times our photometric uncertainty. Furthermore, such a hot white 













dwarf would have a UV continuum that was not detected in Judge et al., who constrain the 
white dwarf’s temperature to less than 10000K if it exists. This low temperature would make the 
magnitude difference even more extreme. The white dwarf possibility is thus definitively excluded. 


m 


On the other hand, the measured flux is consistent with a main-sequence K-type dwarf. Shown 
Figure [3] are spectral models of KO, K5, and K7 stars, with our measured flux shown as a black 
dot. For the 8 And primary, the effective temperature and radius is taken from published results 
(Table [IJ note that the radius is known accurately from interferometry (Piau et ah 2011)). For 
the secondary, the effective temperatures and surface gravities are taken from Castelli & Kurucz 
(2004) and the radii are taken from Cox (2000). While the formal photometric error is smaller 


than the size of the datapoint, lack of knowledge about the companion’s radius and temperature 
make it impossible to give a completely specific spectral classification; the best we can say is 
that the companion is most likely a K-type dwarf. Making a more accurate measurement of the 
secondary spectral type is possible in principle. The simplest way would likely be a similarly 
precise coronagraphic measurement in J band, as the J-K colors of K dwarfs change by about 200 
millimags over the spectral type. Alternatively, an AO-fed integral field spectrograph might be 
able to measure the spectral type from the CO band at approximately 2.3 /jrn. 

The conclusion that the companion is K-type is mostly consistent with previous work. As 
mentioned before, Judge et al. concluded that a main sequence companion would have to be a star 
later than G-type. A K-dwarf has a mass of between 0.6-0.8 M 0 ; taking values of 1.1 - 1.2M 0 for 


8 And A gives 0.3-0.4 as the mass fraction, reasonably consistent with the value of Gontcharov & 


Kiyaeva (2002) of 0.5 ± 0.1. However, there are some discrepancies. If the mass of delta Andromeda 


is so low, the age must be much larger than the reported value of 3.2 Gyr from Decin et al. (2003), 


as sufficient time would not have passed for the star to evolve off the main sequence. This is not 


totally inconsistent, as Decin et al. (2003) express very low confidence in the accuracy of the age 


measurement for this particular star. In light of this, we suggest that the age needs to be revised 
significantly upwards, to over 6 Gyrs. Another more remote possibility consistent with our data is 
that the companion is actually a multiple itself, such as two even lower mass dwarfs. 


The results presented here demonstrate the potential of high contrast imaging applied to 
medium to long-period spectroscopic binaries. In particular, the contrast differences between main 
sequence stars in binaries are readily accessible to a coronagraphic system, and the information 
gained can improve orbit characterization, or as in our case, distinguish quickly between different 
companion possibilities. 
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